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We study the spherically symmetric collapse of a real, minimally coupled, massive scalar field in an 
asymptotically Einstein-de Sitter spacetime background. By means of an eikonal approximation for 
the field and metric functions, we obtain a simple analytical criterion — involving the physical size 
and mass scales (the field's inverse Compton wavelength and the spacetime gravitational mass) of 
the initial matter configuration — for generic (non-time-symmetric) initial data to collapse to a black 
hole. This analytical condition can then be used to place constraints on the initial primordial black 
hole spectrum, by considering spherical density perturbations that re-entered the horizon during an 
early matter-dominated phase that immediately followed infiation. 
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I. INTRODUCTION 

Until now no fundamental elementary spin-0 parti- 
cle has been detected in accelerator experiments, even 
though its existence is predicted by the so far highly suc- 
cessful SU(3)(8)SU(2)(g)U(l) Standard Model 0. Despite 
this lack of knowledge, most models of early universe as- 
troparticle physics include at least one such field, and 
their remarkable agreement with observational data (see 
e.g. 1^ for a recent review) makes pertinent the question 
of whether primordial density fluctuations in the scalar 
field energy density distribution could have collapsed to 
form black holes. 

In particular, several models of inflation predict the 
existence of a period dominated by the energy of a real 
massive scalar fleld, just after the end of the inflation- 
ary epoch 1^-^. During this period, the field can behave 
like dust and density inhomogeneities can undergo non- 
linear growth, which may lead to the formation of primor- 
dial black holes Q. It is during this early matter domi- 
nated phase that black hole formation — via gravitational 
collapse — is the most abundant, and we should therefore 
expect important constraints on inflationary models from 
the overproduction of primordial black holes during this 
epoch. 

In this paper, we address one aspect of this problem by 
obtaining analytical conditions for black hole formation 
during the early matter dominated phase; the constrain- 
ing of the primordial black hole energy density spectrum 
and its implications for inflationary scenarios will be dis- 
cussed in detail elsewhere j?). 

Our model consists of a real, spherically symmetric 
massive scalar field, minimally coupled to gravity, evolv- 
ing in an asymptotically Einstein-de Sitter (EdeS) space- 
time background. The Einstein equations for gravity cou- 
pled to the massive scalar field are solved analytically, 
using an asymptotic expansion for the field (and metric 
fimctions), in terms of the field's Compton wavelength. 



to explore the large mass limit. In such a limit, 
the scalar field behaves like general inhomogeneous dust 
§ , described by the Tolman-Bondi metrics [p|-[l2[ . From 
generic initial data, spherical dust collapse always pro- 
ceeds to a black hol^ and thus, by imposing the condi- 
tion that the eikonal approximation holds, we obtain a 
sufficient criterion for black hole formation from massive 
scalar field collapse. 

This paper is organized as follows: Section II de- 
rives the field equations for the massive Einstein-Klein- 
Gordon system, which are then analytically solved using 
an asymptotic series (in for the field and metric 

functions. The large mass limit of the eikonal approxi- 
mation is the Tolman-Bondi family of metrics, described 
in Sec. III. In Sec. IV, the gravitational collapse of spher- 
ical density perturbations in the EdeS universe — a par- 
ticular case of the general Tolman-Bondi spacetimes — is 
analyzed in terms of initial data. In Sec. V, analytical 
conditions for the initial data to collapse to a black hole 
are obtained by enforcing the validity of the WKB ap- 
proximation to second order in Section VI discusses 
how this condition for black hole formation might affect 
the constraining of the initial (at the early matter dom- 
inated phase) primordial black hole spectrum, and thus 
the inflationary models responsible for generating the ini- 
tial density fluctuation spectrum. Section VII concludes 
with a summary and a brief discussion on avenues for 
future work. 

Geometrized units, in which G = c ~ \, are used 
throughout. 



'Central naked singularities can form from regular ini- 
tial data, but they are not generic: only one non-spacelike 
geodesic can escape from the singularity, which, although 
strong, is also massless (nsl. 
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II. THE MASSIVE EINSTEIN-KLEIN-GORDON 
SYSTEM 

We consider a general spherically symmetric metric, 
written here in normal Gaussian coordinates {t, r, 9, ip}: 

where A = A(i,r), R = R{t,r), and d9? ^ d9'^ + 
sin^ 9d(p^ is the canonical metric of the unit 2-sphere. 

The independent non-vanishing Einstein tensor com- 
ponents are 



Gtt = R-^[-Re^^{2R'A' + 2R" + R-^R'^) 

-2RAR+l+if], 
Grt = -2R-~^{R' + R'k), 



Gag 



-2 



-R 

„:„-2 



e-^^{2RR + R^ + \) - R'' 



=R[RK + R'A 



I A '„2A 



+R!'e^^ - R + AR- K^R], 



(2) 
(3) 

(4) 
(5) 



where the overdot and prime denote partial differentia- 
tion with respect to t and r, respectively. 

For the matter content, we consider a real minimally 
coupled scalar field (j) of mass /i, governed by the Klein- 
Gordon equation: 

ia-^l^)^^o. (6) 

With the spherically symmetric metric (|l|) we have 

4> - e^^ct)" + /i^^ + (A + 2R-^R)<i) 

-e^^{A' - 2R-^R')c^' = 0. (7) 

The stress-energy tensor of the scalar field is given by 

(8) 
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Tab = Va(t>Vb(t> ~ ^gab{Vc(t>V''(f> + fl^cf)'^), 

with the independent non-vanishing components 
1 



Trt — 



-'^ ;2 I 1 2A/,/2 1 2/2 
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Tee - sin-^ ^T^^ = i^R {4> - & <P - M </- ) 
By defining the auxiliary functions 



.2A 



m(t,r) EE -R{R^ + k), 



(9) 
(10) 
(11) 
(12) 

(13) 
(14) 



Einstein's equations can be recast in terms of the first 
derivatives of these two functions: 



k' = -8nRR'{Ttt + T;) - 2R'{R + AR), 
k = 8ttRR'T[, 
m! ^ AttR^R'Tu - AirR^RTrt, 
m = AirR^R'n - AttR^RTI. 



(15) 
(16) 
(17) 
(18) 



Since there are only three independent functions to be 
determined and four equations, only three of these are 
independent, with the remaining one being a constraint. 
Since the scalar wave Eq. (^) is implied by the Einstein 
equations (it follows from the contracted Bianchi identi- 
ties), we shall take it together with Eqs. ( p^ and ( p^ ) 
as our complete set. We take Eq. (^ as the constraint 
equation, since it provides a simple relation between the 
initial data and the initial mass profile. 

To facilitate the resolution of the field equations, we 
introduce a WKB approximation for the field in the large 
mass (/i) limit — when the Compton wavelength of the 
scalar field, is much smaller than the radius A of 

the spherical region where the field is non- vanishing: 



A/i > 1. 



(19) 



In such a limit, we expect wave-like solutions with slowly- 
varying (with respect to t) amplitude: 

(l>{t, r) = /i"^*(^, r) cos ^J,t. (20) 
The stress-energy tensor components are then 

Ttt = - i/i-i$$sin 2/i<-K i($2 +e'*'^$'2) 

x(l + cos2^i), (21) 
Trt = -^M^^'J"!'' sin 2fit + ^fr'^<^<^'{l + cos 2/xi), (22) 
r; = cos 2fit- i^"i$$sin 2nt 
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Tee = ^i?^[-*^cos 2^<-^"^$$sin 2fit 



(23) 



(24) 



The form of the above equations suggests a trigonometric 
expansion for <i>(t, r) of the form 



* = *o + ™(*mn cos nfit + sin niit). 

(25) 



m— 1 n—1 



The other metric functions are expanded analogously. 
This expansion in inverse powers of the mass leads to 
an asymptotic series |lj,|l^ in /i^^, and hence it gives an 
exact solution in the infinite mass limit and approximate 
one otherwise. This method is usually referred to as the 
Lagrange method or the method of averaging ||l^ . 
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Substitution of the expanded metric functions and 
their derivatives in Eqs. (||), jl^), ( p^ and jl^ ) fixes the 
coefficients of ah the trigonometric terms at each order in 
Up to 0(/i~^) we find that the only non- vanishing 
terms are: 

$(t,r) = $o(t,r)+0(Ai-3), (26) 
m{t, r) = mo{t, r) + ii^^mi2{t, r) sin 2iit 

+Ai"^m22(t, r) cos 2^t + 0{p^^), (27) 
k{t, r) = ko{t, r) + /i"^fc22(i, r) cos 2/ii + 0{n~^), (28) 
r) = Ro{t, r) + Ai"^i?22(t, cos 2/ii + CI(m~^). (29) 

Taking the first derivatives of the above equations and 
comparing the leading order terms to the right-hand-side 
of Eqs. (§), (|l|), ^) and (|l|), yields 

mo = 0, (30) 
fco = 0, (31) 
Rl^2m„Ro^ ^k„. (32) 

The leading order terms from Eq. ( [T^ ) give 

Too 27ri?^i?^,$g. (33) 

The first three equations restrict this class of metrics to 
the Tolman-Bondi family, which describes the collapse 
of general inhomogeneous dust. We have thence con- 
cluded that sufficiently — in a sense to be defined precisely 
below — massive scalar fields behave like a dust. 

Spherical dust collapse is pressureless and thus always 
proceeds to a black hole. Therefore, by guaranteeing that 
the field behaves like dust until the complete formation 
of the event horizon, we can obtain a sufficient condition 
for black hole formation from massive scalar field col- 
lapse. This criterion is explicitly obtained by examining 
the next order terms in the WKB approximation. From 



Eqs. (|l|), (^ and (|T|), we have 

fc22(t, r) = 2T:R'o{Ror\l - fco)$o$o, (34) 

mi2{t,r) ^ttRIRo^I (35) 

i?22(i,r) =-|i?o$o- (36) 



These terms have the convenient feature of being simple 
algebraic functions of the known zeroth-order terms and 
their derivatives. The asymptotic expansion guarantees 
that the WKB approximation will remain valid whilst 
all the next order terms remain small compared to the 
leading order ones: 

-2^22 . 1 , 2I-R22I , 1 , -1^12 1 

M 1— < and < - and ^ < -. (37 

fco 2 Rf) 2 Too 2 

The value 1/2 was chosen arbitrarily (it has to be smaller 
than unity!) and changing this value amounts to a change 
in fi. The WKB approximation will therefore breakdown 
when one of the above inequalities is no longer satisfied. 



The asymptotic expansion implies that the largest cor- 
rection will come from the TO12 term, and we have verified 
this numerically. Thus, the validity of the WKB approx- 
imation is given by the condition [where Eqs. (|lj) and 
( p5| ) were used]: 

RorriQ < AiTOoi?o: (38) 

which defines a region on the t — r (hence t — R) plane, 
as shown in Fig.l. The WKB approximation will hold at 
all points in whose causal past the inequality is always 
satisfied. On the t—R plane, it is the set of events outside 
the curve that saturates inequality (p8|), whose past null 
cone is tangent to it. 




WKB approx. 
breakdown 



2M 2M 



FIG. 1. Spherically symmetric dust collapse in {t, R} coor- 
dinates. The solid line is the sphere's surface, and the dashed 
curve that asymptotes R = 2M is the event horizon. The 
patterned region inside the event horizon is where the WKB 
condition breaks down. The eikonal approximation holds for 
events A and B, but not for event C: its past light cone in- 
tersects the region where the WKB condition breaks down. 



III. THE TOLMAN-BONDI LIMIT 

From Eqs. (|l|), (|l|), (|l|), (|l|), we have (hereafter 
dropping the 'O'-indices for simplicity): 

ds^ - -de + ?^^^^dr^ + R\t, r)dn\ (39) 
1 - k{r) 

9 . ^ 1 to' , ^ 

This is the Tolman-Bondi class of solutions. Included in 
this class are the Schwarzschild metric [m{r) — const.], 
the Einstein-de Sitter universe [ R{t,r) = a{t)r oc t^/'^r 
and fc(r) = 0], and the closed Friedmann universe [k{r) oc 
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r ]. The metric is written in comoving coordinates and 
the model describes the pressureless collapse of spherical 
dust shells. Each shell is labelled by the radial coordinate 
r and has a surface area ATTR^{t,r). In the context of 
Tolman-Bondi metrics, shell crossings — the overlapping 
of neighboring shells at finite proper radius — are defined 
by the locus of events for which R > and R' = 0. It 
is clear from Eq. (^ij) that the energy density diverges 
at shell crossing singularities. Even though some curva- 
ture components may diverge at shell crossings the 
spacetime is geodesically complete — analytical continua- 
tions of the metric can always be found (in a distribu- 
tional sense) |Q — and thus they are not real physical 
singularities; rather, they merely signal the intersection 
of matter fiow lines at a given spacelike surface, and thus 
lack of integrability of the field equations beyond that 
surface (it is worth pointing out that they also occur in 
inhomogeneous spherical Newtonian collapse) . 

Assuming R' ^ 0, we can integrate Eq. (^0|) paramet- 
rically to obtain 



tiVir) =^ta (r) + (77 + sin 77) , 
R{T],r) 



2m 



2 V 



(42) 
(43) 



where to{r) is an arbitrary constant of integration (to be 
fixed by the initial data), and < ?/ < tt. The physical 
initial data consists of a mass distribution m(r) and a 
velocity profile R{0,r). There is gauge freedom for the 
scaling of r and we shall fix it by requiring that it coin- 
cides with the initial area radius: 

R{0, r) = r. (44) 

A. Time-symmetric initial data 

Considering time-symmetric initial data, 

7?(0,r) = 0, (45) 

implies, from Eq. (40), fc(r) — 2m /r and tQ{r) — 0. 
Equations (p2|) and (43) then simplify to 



(46) 
(47) 



= ^ (?7-fsin77), 
R{-q, r) —r cos^ ^. 



A shell with initial proper area Aur^ will thus collapse to 
vanishing area radius in a (comoving) time 



(48) 



For inhomogeneous mass distributions {m 7^ const, x r^), 
different shells become singular at different times; in the 




homogeneous case, all the shells collapse to zero area ra- 
dius at the same time [ ]T7| . 

A particular solution is specified by an initial mass 
function m[r), given uniquely by the initial density dis- 
tribution: 

m(r)^2TT $^(0,f)f2df = 47r / Ttt{0,r)f'^df. (49) 
Jq Jo 

If the mass function m(r) approaches a constant value 
AI when r — > +00, then AI is the ADM mass of the 
spacetime. 



IV. THE EINSTEIN-DE SITTER UNIVERSE 

The EdeS universe is a particular member of the 
Friedmann- Robertson- Walker metrics[]. 



a'{t) 



dr^ 



1 - Kr^ 



r'^dVl? 



(50) 



where K = —1, 0, or -1-1, for a hyperbolic, fiat or spherical 
spatial geometry, respectively. The Einstein equations 
are simply 



P = -3(p + p)-, 
a 

d At: , 
a 6 



Stt 



K 
7^' 



(51) 
(52) 

(53) 



where p is the proper energy density and p the pressure. 
For pressureless matter distributions (p = 0) and vanish- 
ing spatial curvature {K — 0) — the EdeS universe — the 
solution to the Friedmann equation (p3|) is 



a{t) ~ aot3 , 



(54) 



where ag is an integration constant. This model de- 
scribes an open geometry (the K — spatial sections 
are diffeomorphic to R^) in the presence of a constant 
non-zero energy density distribution. Even though the 
EdeS spacetime is conformally fiat, its causal structure 
is quite different from asymptotically fiat geometries. In 
particular, and unlike Minkowski or Schwarzschild, past 
null infinity for EdeS is spacelike, and thus past particle 
horizons exist [which can be observer (i.e., r) dependent]. 

The EdeS spacetime is a good approximation to the 
large scale structure of the universe during a matter 
dominated phase, when the averaged (over space and 



^ Which are a subclass of the Tolman-Bondi metrics [cf. Eq. 
(Ei)], with R{t,r) = a{t)r and k{r) = Kr^. 
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time) energy density evolves adiabatically and pressures 
are vanishingly small, as, e.g., immediately after infla- 
tion [§1). We shall therefore adopt such a metric to 
model the collapse of overdensity perturbations in the 
early matter dominated phase that followed inflation. 




FIG. 2. Conformal diagram for the EdeS spacetime. The 
dashed vertical line on the left denotes the coordinate (cen- 
tral) singularity. PNC is the past null cone of a timelike ob- 
server O at event p. Worldlines starting from (spacelike) null 
infinity outside p' s PNC are causally disconnected from O at 



A. Dust collapse in the Einstein-de Sitter universe 

Our model consists of a general Tolman-Bondi metric 
and a spherically symmetric matter distribution of finite 
proper radius -R(0, r) = r = A, with an overdensity (.{r) 
with respect to the constant background energy density, 

Pc- 



T„(0,r) = p(r) = pe[l + CM], 



(55) 



with i^(r) — ^(r)O(r — A), where C{r) is an arbitrary real- 
valued function and Q is the unit Hcaviside function. For 
r ^ A, the metric asymptotically approaches the EdeS 
metric and the finite background energy density pc im- 
plies the existence of a cosmological horizon with proper 
radius = (87rpc/3)~^/^. Clearly, the existence of 

such a horizon (or, equivalently, expanding exterior ge- 
ometry) precludes the use of time-symmetric initial data. 

Let us then define the generalized Hubble parameter 
by 



iJ(i,r)^|. 



(56) 



Then, from Eqs. (HC 



we have 



2m % 
r" 2 ' 
2m 2 '^o 



^oM-^Bin^:^ 

k{r) 



— cos „ 
r 2 



2 m 

r 



1 - 



2m 



to{r) 



m 
p72 



[?7o(r) -f-sin%(7-)]. 



(57) 
(58) 
(59) 



where Ho{r) = H{Q,r) — R{0,r)/r. Initial data consists 
of a mass profile m(r) and an initial Hubble parameter 
Ho{r). We shall consider here the case of an EdeS uni- 
verse with initial expansion rate 



Stt 



(60) 



From Eqs. ( }49| ) and (p5D, the mass function is then 
m(r) = 47r / Tu{0,x)x^dx = —py{l + ^) 



M 





r\3 



(61) 



where ^ is the volume average of ^(r) and M is the mass 
inside a sphere of proper area 47ri?^(i, A). 




Surface of density perturbation, S 



Future null cone of S 



FIG. 3. Conformal diagram for spherical dust collapse in 
the EdeS universe. As the overdensity collapses (dashed 
curve), it first crosses the future cosmological horizon (FCH) 
of an r = observer, and, at a later time, the black hole 
event horizon, eventually ending up at the R = central sin- 
gularity. Events outside the FCH will never be seen by the 
r = observer. Events outside the future null cone of S (the 
overdensity's spacelike 2-surface) will never be influenced by 
the overdensity perturbation. 



V. CONDITIONS FOR BLACK HOLE 
FORMATION 

We want to impose the condition that the stationary 
phase approximation holds [cf. Eq. (Ssl)]: 



Rm' < p,R'm, 



(62) 



where 
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R' 



tan 



^ ~k 
R = —Vkt&n ^, 

to' 3 fc' 
^(") = TO - 2fc- 



(64) 
(65) 
(66) 



For simplicity, we consider the energy density distribu- 
tion inside the overdense region to be constant and hence, 
from Eqs. (^, ( |65| ) and (66), we have 7(r) = io(f) = 0. 
The relevant WKB condition (p3) becomes then 



2M 



< fi cos^ - cot 



where 



'1 



G[0,1) 



(67) 



(68) 



is a monotonically increasing function of the magnitude 
of the volume-averaged overdensity, satisfying p(0) = 
and lim^^_|_3^ P = 1- 

Outside the overdensity region, for r > A, the Tolman- 
Bondi solution is exact: the EdeS metric is a particular 
case, trivially obtained from the general Tolman-Bondi 
form by setting R{t,r) — {t/tCj'^/^r and k{r) — 0, where 
ti is the initial slice. Therefore it suffices to enforce the 
above inequality (|6^) for < r < A, to ensure that the 
WKB approximation will hold for all r > 0. If all the 
shells with < r < A cross the event horizon before the 
WKB breaks down, then we can conclude that a black 
hole has formed. 

In order to find the loci of events defined by the WKB 
condition and place them wholly inside the event horizon 
(EH) , we need to (i) find a suitable parameterization for 
the curve that bounds the region defined by Eq. (67), 
and (ii) solve for the EH with the same parametrization 
for the null generators. 

The natural such parameterization for the curve that 
saturates inequality ( |67|) is 77* (r) [ since rj and p (through 
^) depend solely on r], given by: 



2M 



2 



^P"\l^ = AfCOS" y cot 2 



(69) 



From Eqs. ( |39| ) and ([42|), we can parameterize the 
null generators by ?7eh, obtained by integration of the 
outgoing radial null geodesies equation: 



dr 



{i?'(fc-i - 1)-^ - [t'oVk + ^7(77 + sin?;)] } R.-\ 

(70) 



For r < A, this simplifies to 
d?7EH ^ A^ 



dr 



(71) 



which is trivially integrated to give 
VEnir < A) = ?;eh(0) + sin"-^p 



2M r 



(72) 



Outside the overdense region, 7^0 and ^ 0, which 
makes it impossible to integrate the radial null geodesies 
equation analytically. Since the exterior geometry is not 
asymptotically flat, we cannot use Birkhoff's theorem to 
place the EH at i? = 2M. We can, however, solve for the 
apparent horizon (AH) — the outer boundary of a closed 
spacelike surface whose normal null geodesies have van- 
ishing divergence — at any given t = const, hypersurface. 
In spherical symmetry, the equation for the AH is simply 
the requirement that a surface R(t, r) = const, becomes 
null: 

R'^R.a ^ 0. (73) 
From Eqs. (^) and (^o[), we have then 

i?(tAH,r) = 2TO(r), (74) 
which, via Eq. (64), defines a curve ?7AH(f'), given by 



„ -1 r 2M 
2 cos- -^—p. 



(75) 



Since the AH must be contained or coincide with the 
(spacelike section of the) EH, we must have ?7ah(A) > 
?7eh(A), which gives an upper bound for r7EH(0) = t/q: 



2M 



770 < 2 cos \ -—p - sin \ -—p. (76) 



A 



2M 



A 



This upper bound on 770 guarantees that the approximate 
EH thus estimated coincides or lies inside the actual EH; 
hence, if the WKB approximation holds inside the for- 
mer, it will obviously hold inside the latter. 

We now need to guarantee that the breakdown do- 
main of the WKB approximation is placed entirely in- 
side the EH. From Eq. (|69|), 77* is implicitly defined by 

an expression of the form /[77*] = ^P^^\J^^> where / 
is a monotonically increasing functional of 77, ; therefore, 
V* > ??eh(A) /[77*] > /[7;eh(A)]. This enables a suffi- 
cient condition for the WKB approximation to hold (in- 
side the true EH) , by imposing a lower bound on rj^ : 



V* > '7eh(A) > 2cos ^ 



2M 



A 



-P- 



(77) 



From the above equation, together with Eq. (|67|), we 
finally obtain 



/iM > 



1 - —p^. 



(78) 



Matter configurations with parameters satisfying the 
above inequality will collapse to form a black hole. This 
is a sufficient, but not necessary, criterion, due to both 
the adiabatic ansatz {fiX 3> 1) and the approximation for 
the EH. There are two limits that should be pointed out: 
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A.l. p -> 1 



B.2. HoX ~ 1 



In the hmit of very large overdensities, the inequal- 
ity approaches the one for the asymptotically flat case 
1^, as expected. The simple physical reason is that the 
ratio Pc/Pport. <C 1, and therefore having a vanishing 
or finite background energy density becomes irrelevant; 
the overdense region does not "feel" the exterior geome- 
try. This can be heuristically understood as follows: The 

1/2 

typical cosmological expansion rate is {R/R)cxp ^ Re , 
whereas the collapse rate is {R/R)co\\ ^ ~[(1 + 0/Oc]^^^^- 

Then, ^ = |(i),,n|/(^)e.p ^ + » 
expansion that the overdensity would have experienced 
during a (proper) time fcoU, had it not collapsed during 
that time, is negligible compared to its initial physical 
radius. 



If ^ > {HqX)^^, then the right-hand-side of the 
inequality becomes imaginary and the approximation 
breaks down. Physically, this corresponds to an overden- 
sity that becomes so large that a spacelike hypersurface 
would curve upon itself and form a disjoint spatial uni- 
verse, which does not correspond to a black hole [^ . 
Provided ^ < {HqX)^^ < 1, this limit is equivalent to 
the p -^1 limit previously described, and thus still valid. 
From an astrophysical viewpoint, this is the most inter- 
esting limit, as it constrains the density perturbations at 
horizon crossing. 



VI. CONSTRAINING THE PRIMORDIAL 
BLACK HOLE SPECTRUM 



A. 2. p^O 

When the overdensity vanishes, although the right- 
hand-side of the inequality approaches 3/2, this result is 
meaningless because it was based on an estimate for the 
EH location that becomes totally unreliable: ?7eh(0) = tt, 
when p = 0. This is trivial, since p — corresponds to an 
unperturbed EdeS universe. For finite p ^ 1, provided 
the spatial extent of the overdensity is sufficiently large, 
A > 2Mp^, the inequality is still valid. 

We can also look at this condition in terms of the cos- 
mological horizon, Hq^. Noting that 2M = H^X^{1 - 
p'^)^^ and recalling that p = {1 + ^)~^^'^, we can rewrite 
the condition for black hole formation as 

^,M > yi-{HoXn. (79) 

The two relevant limits are now: 

B.l. HoX<^ 1 

The length scale A is much smaller than the horizon 
radius Hq^, and f can be sufficiently large to recover the 
p — > 1 limit and thus the asymptotically flat results. This 
limit also provides a useful self-consistency check for the 
validity of the WKB approximation. Let us take HqX ^ 
1 and f ~ 0(1); then, the inequality becomes H^X'^fi > 
|. Since, by hypothesis, HqX ^ 1, the necessary and 
sufhcient condition to enforce the inequality is pA 3> 1, 
which is precisely the adiabatic condition on which the 
whole argument is based. Thus, we confirm that if the 
WKB approximation is valid initially, it will remain valid 
up to the formation of the EH. 



Even though inflation is still a scenario and not 
yet a paradigm for cosmology, its basic predictions 
appear fairly robust in that not only they solve the 
("naturalness" and "fine-tuning") shortcomings of the 
"old" standard Big Bang model, but they also present 
a viable physical mechanism for the primordial density 
fluctuations that seeded the subsequent structure forma- 
tion. 

As previously mentioned, several inflationary models 
predict the existence of a period dominated by the energy 
density of a real massive scalar field, just after the end 
of inflation, hereafter denoted by tend- During this stage, 
the field undergoes coherent oscillations, behaving effec- 
tively like a dust, from tend until ^rh ("reheating"), when 
it rapidly decays into relativistic particles . Any den- 
sity perturbations crossing the horizon during this dust- 
like phase will inevitably (there is no pressure) collapse to 
form black holes, provided they are sufficiently spherical, 
as aspherical growth precludes the gravitational collapse 
of the inhomogeneities j22j. Clearly, it is during this in- 
termediate matter-dominated phase that black hole for- 
mation via gravitational collapse of (spherical) overden- 
sities is the most abundant. We should therefore expect 
important constraints on inflationary scenarios (e.g., on 
the decay width of the scalar fleld, F) from the over- 
production of primordial black holes (PBH) during this 
period. 

The present process will only be relevant for black hole 
formation if the time interval, F~^, after which the field 
decays into radiation is larger than the collapse time 
to the central singularity. We note that it is sufficient 
to require that the scalar field dominates just until the 
overdense region has collapsed through its Schwarzschild 
radius, since a subsequent decay into radiation (before 
vanishing proper area is reached) would not destroy the 
already formed black hole; furthermore, the proper time 
to reach the formation of black hole horizons (i.e., the 
proper time elapsed between horizon crossing and the 
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emission of the null generators) is much larger than the 
time it would take for a black-hole-sized configuration to 
collapse to a singularity. 

Since spherical symmetry is assumed throughout, the 
necessary and sufficient conditions for black hole forma- 
tion from a massive scalar field configuration with pa- 
rameters {/i, M, A} are that black holes: 



have time to form: iuc + ^coU < ^rh, 



(80) 



will form: /iM > -^J I - (i/oA)^^, 



31^ 



where in the first inequality, the subscripts refer to "hori- 
zon crossing", "collapse" and "reheating", respectively. 
Because the synchronous gauge is being used, the three 
timescales can be trivially compared, thus justifying in- 
equality (|^. Since we are considering homogeneous den- 
sity distributions, the proper time for the collapse of the 
overdense region will depend solely on its density. Setting 
7? = TT, from Eqs. (B|), (E|) and (|57l)-(|9|), we obtain 



icon - ii/o"^r^(l + 0{tt - sin^i x - x), (82) 



(83) 



For large overdensities, f ^ 1, ^coii ^ k'^-^o ^^^i 
whereas in the limit of small perturbations, ^ <C 1, the 
collapse time is 



(84) 



In the EdeS universe, Hq = ^ ^ y/Snpc/S. At 



horizon crossing, HX — 1; thus 

2 2 _i 
Hc 3 3 



(85) 



With Eq. (|84|), we can rewrite Eq. ( pOD as (dropping the 
bar for simplicity) 



(1) 



(86) 



That is, at a given slice t = ^Hq^, there will be a mini- 
mum fractional density perturbation at horizon crossing, 
below which the scalar field will decay into radiation be- 
fore the overdensity has time to undergo complete grav- 
itational collapse. 

Let us now examine the WKB condition at horizon 
crossing. The mass of a perturbation with fractional 
overdensity at horizon crossing, is 

Af = ^7rpc(l + OA' = ^^o(l + OA' 



(87) 



Equation (pl[) becomes then 



which places another lower limit on ^: 

e > 1 - If^'M' ^ 



(88) 



(89) 



We have then two independent lower bounds on ^, that 
tell us, respectively, that black holes can and will form. 
If > |, then the dust approximation remains reliable 
and Eq. (^) is always satisfied. In this case, Eq. (|8^) be- 
comes the relevant constraint, and the crucial quantity is 
Hat^H, which, although model dependent (through tpiH, 
or, equivalently, T"^, for a given horizon crossing time), 
is expected to be large, i?oiRH > 10® [^ . 

If fj,M < 3/2, then the dust approximation may be- 
come unreliable and the fractional density perturbations 
are of order unity, therefore departing from the linear 
regime and being heavily suppressed in the initial per- 
turbation spectrum, in which case the dust-like phase is 
irrelevant for the production of PBH from the gravita- 
tional collapse of massive scalar fields. 

We shall then consider only configurations satisfying 
fiM > 3/2. Let us first rewrite the density contrast in 
more familiar notation: 



Pc Pc 



(90) 



where 6 is to be evaluated at horizon crossing. 

For the spectrum of density fluctuations crossing the 
horizon, we consider a Gaussian probability distribution 
(as predicted by inflation pH): 



PiS) = 



27ra 



cxp 



"2a2 



(91) 



where a is the root-mean-square amplitude of the density 
fluctuations. Provided we restrict ourselves to linear per- 
turbations, the Fourier modes will evolve independently, 
thereby preserving Gaussianity. (Non-Gaussian effects 
arise if (5 3> cr, when linearity is spoiled by mode-mode 
coupling, but their influence on the initial PBH mass 
function for a dust-dominated era has been shown to be 
neghgible ||). 

The quantity of interest in the constraining of PBH 
production is the mass fraction of the universe that col- 
lapsed to a black hole during the intermediate dust-like 
epoch: 



(92) 



where ppbh is the initial (i.e., formed during that epoch) 
energy density of primordial black holes, and pc is the 
background energy density, ^min is given by Eq. ( |8^ ) 
and the upper integration limit is imposed by linearity. 



8 



Expressed in this way, (3 is model dependent, since ^min 
depends on Hq and tRH (equivalently, r~^). 

(3 can be observationally constrained via (i) the re- 
quirement that the present density fraction of PBH's, 
^^PBH < 1; and (ii) hmits from the Hawking radiation, 
from upper bounds on the PBH population at subsequent 
epochs, e.g., baryogenesis and nucleosynthesis. These 
constraints on (3 will lead to constraints on the cosmolog- 
ical parameters (3 is implicitly dependent on, namely the 
decay width of the scalar field and the time for horizon 
crossing. Work in this direction is currently underway 
and will be presented in detail elsewhere H . 



VII. CONCLUSIONS 

We have shown that there is an adiabatic regime 
(A/i 3> 1), in which massive scalar fields in an asymp- 
totically EdeS background behave like a dust and there- 
fore inevitably collapse to form black holes. Enforcing 
the adiabacity condition to second order in lead to 
a simple analytical criterion to be obeyed by the initial 
data to collapse to black holes. This criterion is valid 
for spherical density perturbations of any finite physical 
radius A ; in particular, they can equal the horizon size, 
H\ — 1. Perturbations with spatial radius larger than 
the horizon radius are physically uninteresting because 

(i) the fractional density perturbation required to sat- 
isfy the WKB condition is very small, ^ <C (HX)^^, and 

(ii) they are causally disconnected from the observable 
universe. 

Immediately after inflation, there can exist a short pe- 
riod dominated by the energy density of a real massive 
scalar field |^-^. By considering spherical density per- 
turbations that re-enter the horizon during this epoch, we 
obtained a lower limit on the fractional density pertur- 
bation (at horizon crossing) that would allow black holes 
to form before the field decays into relativistic particles. 
This lower limit can then be used to compute the ini- 
tial mass function of PBH at the early matter-dominated 
phase. Constraints on the PBH mass spectrum from its 
subsequent evolution will ultimately lead to constraints 
on some free parameters of particle astrophysics models 
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